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Abstract 

Network virtualization techniques allow for the coexistence of many 
virtual networks (VNs) jointly sharing the resources of an underlying sub¬ 
strate network. The Virtual Network Embedding problem (VNE) arises 
when looking for the most profitable set of VNs to embed onto the sub¬ 
strate. In this paper, we address the offline version of the problem. We 
propose a Mixed-Integer Linear Programming formulation to solve it to 
optimality which accounts for acceptance and rejection of virtual network 
requests, allowing for both splittable and unsplittable (single path) rout¬ 
ing schemes. Our formulation also considers a Rent-at-Bulk (RaB) model 
for the rental of substrate capacities where economies of scale apply. To 
better emphasize the importance of RaB, we also compare our method 
to a baseline one which only takes RaB into account a posteriori, once 
a solution to VNE, oblivious to RaB, has been found. Computational 
experiments show the viability of our approach, stressing the relevance of 
addressing RaB directly with an exact formulation. 

Index Terms Virtual Network Embedding, Network design, Mixed- 
Integer Linear Programming 

1 Introduction 

According to much of the recent literature, network virtualization techniques 
are becoming one of the distinctive features of the new generation of network 
architectures [T]. In nuce, network virtualization consists in decoupling the 
traditional role of current generation Internet service providers into two inde¬ 
pendent roles: management of the physical network and service provisioning. 
Two new actors are thus identified: Service Providers (SPs), who aggregate 
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physical resources so to realize a set of Virtual Networks (VNs) by which end- 
user services are provided, and Infrastructure Providers (InPs), by whom the 
physical network is managed and from whom physical resources are rented. For 
a survey on the topic, we refer the reader to [5]. 



Figure 1: An embedding of two VNs onto the physical substrate. 

The Virtual Network Embedding problem (VNE) arises so to efficiently allo¬ 
cate the resources of the physical layer to the various VNs. It calls for a mapping 
of nodes and links of (a selection of) virtual networks onto nodes and links of the 
physical layer (see Figure]^ for an illustration) so to maximize a profit function. 
From an optimization point of view, the key aspects of VNE are the following: 

• online VS offline: depending on the situation that we consider, VN re¬ 
quests can either arrive dynamically over time or be known a priori^ 

• resource requirements, holding for each VN request (usually of node CPU 
power and bandwidth); 

• locality requirements, imposing that a virtual node be mapped only to a 
specified subset of the physical nodes, (e.g., those belonging to a given 
region); 

• admission control, i.e., the possibility of accepting or denying a VN request 
(e.g., because it is not sufficiently profitable or too resource demanding); 

• routing schemes, depending on the architecture and protocols (e.g., split- 
table routing for UDP traffic or unsplittable single path routing for MPLS). 

Eor both the online and offline cases, with and without admission control, 
VNE is A/'U-hard even with splittable routing by reduction from the multiway 
separator problem [5]. VNE is still A/'U-hard even if we ignore the bandwidth 
requirements due to admitting the generalized assignment problem as a special 
case. If admission control has been carried out and the node mapping is given, 
VNE is polynomially solvable with splittable routing, but it is still A/U-hard in 
the unsplittable case, becoming equivalent to the unsplittable multicommodity 
flow problem. 

Most of the literature on VNE adopts a two-phase heuristic, carrying out 
node and link in sequence [HIHIIS]. To our knowledge, the only “partially exact” 
approach is in [7], where a Mixed-Integer Linear Programming (MILP) formu¬ 
lation is proposed to solve a single step of the online problem to optimality, i.e.. 
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to embed a single VN request or a small group of them. The same formulation 
is then used, essentially unmodified in its basic structure, for variations of the 
original problem such as in [8]. 

To our knowledge, the offline VNE problem has been almost disregarded in 
the literature where, in the few instances in which it is mentioned, it is considered 
only as a mean of comparison for the competitiveness of online algorithms. In 
spite of this, we do believe that the offline problem is of large relevance on its 
own. This is the case, for instance, when the VN requests ask for services (e.g., 
online gaming networks), which, after provisioning starts, are supposed to last 
for a very long long time span. Such requests can be issued long before being 
actually instantiated, thus giving ample time for offline planning. 

Differently from previous work, we also speculate that, in the interaction 
between SPs and InPs for the rental of physical resources, economies of scale 
would typically apply. A new aspect of VNE thus arises, with a direct impact 
on the profits which are maximized in the problem: 

• rent-at-bulk resources (RaB): physical resources are rented in bulks, with 
volume discounts. 

In this letter, we address the offline version of VNE. We propose a MILP 
formulation which entails the admission control aspect, allowing for both un- 
splittable and splittable routing schemes, with a RaB scheme for the rental of 
physical resources. Computational experiments are carried out on a dataset 
composed of long-haul and datacenter networks so to assess the viability of our 
approach. We evaluate both the impact of the two different routing schemes 
and as well as of the new RaB aspect. To better illustrate the relevance of the 
latter, we also compare our solutions to a method which, first, solves a version 
of VNE which is oblivious to the costs of capacity installation and, only then, 
computes (a posteriori) the corresponding RaB costs. 

2 MILP formulation 

Let R be the set of VN requests. Let = (V°,A°) be a directed graph 
representing the physical substrate, with total node and link capacities of Bi 
for all i G and Kij for all {i,j) € A°. Let U and Q be the set of capacity 
bulks of different size for nodes and links. Let a“ and for u S C/ and q G Q, 
be the RaB costs for nodes and links. Economies of scales dictate 

Ui — Uj 

oqi OQi 

and ^ ^ for all Ui,Uj G U with Ui < uj and qi,qj G Q with qi < qj, i.e., 

decreasing unit costs for larger bulks of rented capacity. For each request r G R, 
let V’’ be the set of virtual nodes, with node requirements t’’ for each i G . 
Virtual arcs are implicitly represented via the traffic matrix O'" 
where each component is a demand between two virtual nodes v,w G V’’. 
For each r G R and u G V, let V°(r, u) denote the set of physical nodes on 
which the virtual node v can be mapped (due to locality restrictions). For each 
r G R, let p*" > 0 be the corresponding proht. 
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Let the variable j/’' G {0,1} take value 1 if the request of index r G R is 
accepted and 0 otherwise. Let G {0,1} be equal to 1 if the virtual node 
V G oi request r G R is mapped onto the physical node i G V^{r,v) and to 
0 otherwise. Assuming, for now, an unsplittable routing, let take value 1 

if traffic between the two virtual nodes v,w G for request r G R, is routed 
over the arc {i,j) G and 0 otherwise. This way, we have a multicommodity 
flow with a commodity per triple r, u, w. Let the integer variables gf, hfj denote 
the amount of bulk of capacity u G U and q G Q rented on the physical node i 
and physical link respectively. Our MILP formulation reads: 


\ ^ r r \ ' u \ ^ u 

max 2 ^P y - 2^ Qi - 

E/^^ E hi 

(1) 

rGR ueU i£V° 

qeQ {i,j)eAO 


, \ ^ r r 

s.t. 2_^ x^i = y 

{r,v) 

WgR,vGV" 

(2) 

\ ^ \ ^ iT r ^ \ ^ li 

/ ^ ^ ^v^vi — / ^ ^ 9i 

rGRv^V^ u^U 

Vi G 

(3) 

E E 

r^Rv^w^V^ q€Q 

hi V(i,i)GA° 

(4) 

E ^9i < Bi 

u^U 

Vi G V° 

(5) 

qeQ 

V(i,i) G A° 

(6) 

E /s”™- E i 

'r,vw r r 

ji ^vi ^wi 


WgR,v,wG PCi G 

(7) 

r r rr.vw _ c r\ i 1 u i a 

y G {0,1}, 9i ,hjj 

G Z+. 

(8) 


Constraints ([^ enforce that each virtual node be mapped onto a single physical 
node meeting its locality requirements if the corresponding request is accepted, 
and to none otherwise. Constraints ^ and guarantee that the physical 
node and link capacity that is used does not exceed the rented one, while Con¬ 
straints (§ and <§ impose that no more than the total available physical ca¬ 
pacity be rented. Constraints 0 are multicommodity flow balance constraints 
which, differently from those of a standard multicommodity flow problem, have 
a variable right-hand side. This way, physical node i acts as source node if 
xjjj = 1 and x'^^ = 0, as sink node if = 0 and xj^j = 1, and as a regular inter¬ 
mediate node if x^j = = 0. If xJJ, = xj„j = 1, then the two virtual nodes v, w 

are mapped onto the same physical node (the so-called co-location) and, hence, 
their traffic demand does not need to be routed. Constraints Q denote the 
nature of the variables. Due to variable being defined as binary, a single 

path routing is enforced. A splittable routing can be obtained by relaxing 
to be in [0,1] which, this way, denotes the fraction of flow of request r G R 
between v,w G V'' that is allocated on the arc {i,j). 

Note that, although the formulation proposed in [7] can be used in the offline 
setting as well, it does not allow for admission control, it lacks the RaB scheme 
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for physical resources, and it only allows for a splittable routing. Since such 
formulation aggregates, for each VN, all the flow between pairs of virtual nodes 
that are mapped on the same pair of physical nodes, by imposing integrality on 
the corresponding flow variables we would introduce extra routing constraints, 
thus incorrectly forcing all such heterogeneous flows to share the same physical 
path. 


3 Computational comparison 

In this section, we report on a set of computational experiments carried out to 
assess the effectiveness of our MILP formulation when solved via a commercial 
branch-and-bound code, as well as to evaluate the impact of the different aspects 
of the problem on its solvability and on the quality of its solutions. We adopt 
the state-of-the-art MILP solver CPLEX 12.6 with default parameters, halting 
the execution of CPLEX as soon as a solution with an optimality gap < 1% is 
found. The experiments are run on a default desktop computer with an Intel i7 
processor and 16 GB of RAM, within a time limit of 3600 seconds. 

We consider two types of topologies for the physical networks. We take 
five SNDlib [5] instances (abilene, atlanta, france, germainySO, nobel-eu) 
to model “long-haul networks” representing large-scale backbone networks with 
geographically scattered data centers, each represented as a single aggregate 
node, and five transition-stub instances with |E°| and |A°| in (13,14, 23, 31,45), 
(30,48,60,96,148), generated along the procedure in [10], representing “data¬ 
center networks”, that is, clusters of computers connected via short-haul links. 
For each undirected link in the original instance, two antiparallel arcs are in¬ 
troduced (with possibly different capacities), thus creating directed networks. 
Virtual node and link capacities are randomly generated with values equal to 
5, 10, 50, 500, chosen with probability 0.1, 0.4, 0.4, 0.1, respectively. From 
each topology, we generate two different physical network instances by varying 
the random seed, thus creating a total of 20 physical network instances. We 
consider bulks of size 1, 10, 100, with costs 1, 5, 25, for both nodes and links. 

We couple each of these 20 substrates with a set of 10, 15, 20, 25 virtual 
network requests (Req). Each VN has a random number of nodes between 2 
and 10, with a traffic demand with a density of 50%. Node and link demand 
requirements are randomly sampled as in the physical layer, with a scaling of 
0.3, 0.4, 0.5 {Scat). Embedding profits are set to 500. The resulting data set is 
composed of 240 instances, 10 per configuration of the Req and Seal parameters 
(thus representing problems with different load levels) and type of substrate. As 
to the locality aspect, each set V°(r, v) is constructed by first sampling uniformly 
at random a cardinality factor 7 ” from the interval [|,l]- Afterwards, every 
node u € V° is added to V°(r, u) with probability 7 ”, uniformly at random. 

The computational results for VNE with RaB are summarized in the first 
half of Table for both splittable and unsplittable routing schemes. To show 
the impact of RaB on the solvability of the problem, in the second half of 
the table we relax the integrality constraints on the variables and /iT, thus 
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simulating the case of a standard (constant) price for embedding, proportional 
to resource consumption. Note that, this way, the cost of capacity installation 
is underestimated as any optimal solution will only include fractional amounts 
of the cheapest bulk. Further, to stress that, in an RaB scenario, RaB cost 
are relevant and should be accounted for directly in the model, we introduce a 
baseline method (a heuristic for the VNE problem with RaB) which, based on 
the solution without RaB, computes only a-posteriori the corresponding RaB 
costs for capacity installation. 

Each row reports data averaged (in arithmetic mean) over the 10 instances 
with the same parameter values Req and Seal. The total profit for the embedded 
VNs in the best feasible solution that is found (averaged over the 10 instances 
per row) is reported in the Profit column. The columns ff and Time show the 
number of instances that are solved to optimality within the time limit and the 
corresponding average computing time. The Gap column reports the average 
integrality gap for the instances that are not solved within the time limit (or 0 
if all of them are solved). The percentage difference between the exact and the 
baseline methods is reported in the Impr column. Average data are reported in 
boldface. 

The comparison between our exact approach and the baseline method indi¬ 
cates that the former can lead to much more profitable solutions than the latter, 
with an average improvement of 23.4% which, after a closer inspection, can be 
found as high as 50% on 12.5% of the cases. Note that the data center instances 
are harder to solve than the long-haul network ones, yielding, on average, larger 
computing times, larger gaps, and fewer optimal solutions. In spite of this, we 
still obtain reasonably small gaps (11.5% on average) and, most interestingly, 
an even larger improvement w.r.t. the baseline, equal to 27.3% on average for 
the two routing schemes. 

We now focus on the differences between the two routing schemes over both 
types of instances. We remark that, with splittable routing, we can achieve 
larger profits since the corresponding formulation is a partial relaxation of the 
formulation for the unsplittable case (obtained by relaxing the integrality of the 
flow variables). In spite of this observation, the table shows that the profits 
between the two schemes differ, on average, by no more than 2% (both with or 
without RaB), thus indicating that, in practice, the option of a splittable flow 
does not provide any advantage. What is more, for splittable routing, we register 
an increase (w.r.t. the unsplittable case) in the average computing time (for the 
instances that are solved to optimality within the time limit) of a factor of 1.07 
with RaB and of 1.2 without it, as well as an increase in the number of unsolved 
instances by 16 with RaB, going up to 23 without it. This outcome is, most 
likely, a consequence of the splittable case yielding a less structured problem, 
for which finding good quality feasible solutions via the primal heuristics called 
by the MILP solver is not as easy as for the single path case, where all the 
variables are binary, as well as due to the more effective presolve phase that the 
solver can carry out for pure binary problems. 

When removing RaB, we observe larger gaps for the unsolved instances, with 
an increase, for splittable and unsplittable routing, of 3.6 and 1.3 percentage 
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Number of VN requests 


Number of VN requests 


Figure 2: Average profits plotted against the number of requests, with a scaling 
of 0.3, 0.4, 0.5, with splittable and unsplittable routing schemes and with and 
without RaB. 


points. Interestingly though, we observe, on average, a substantial decrease in 
computing time of 53% and 59%, for the two schemes, as well as an increase in 
the average number of instances solved to optimality by, respectively, 2.9 and 
3.2. This is not surprising since RaB introduces a network design aspect, making 
VNE more similar to classical network design problems which are typically hard 
to solve. Nevertheless, as shown via the comparison with the baseline method, 
such aspect should not be neglected. 

To conclude, in Figure we compare the profits for the embedded VNs 
requests in the four configurations (splittable or unsplittable routing, with or 
without RaB), for different values of the scaling parameter Seal (thus progres¬ 
sively increasing the load level of the instances). Interestingly, the charts report 
the same qualitative growth for all the four cases. Nevertheless, we observe that 
the profits for the case without RaB are larger than those obtained for RaB, 
a feature which better underlines how, for a practical scenario with RaB, the 
model without it underestimates, often substantially, the actual costs. 

4 Conclusions and further research 

We have presented an exact MILP formulation for the offline VNE problem, 
calling for the selection of the most profitable subset of VNs to embed onto the 
substrate. Differently from previous work, our formulation is suitable for both 
splittable and single path routing schemes and introduces a new aspect of VNE, 
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the RaB scheme for the rental of capacities from the substrate. Computational 
experiments have shown the importance of RaB when comparing to a method 
which neglect it, as well as the overall viability of our exact MILP approach. 
Future work includes the study of the polyhedral structure of the problem, so 
as to introduce tighter constraints yielding better bounds and, overall, a faster 
solution process. 
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Table 1: Results for splittable and single path routing schemes, with and without rent-at-bulk. 
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